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ABSTRACT 

The origin of the soft X-ray excess emission observed in many type-1 AGN has 
been an unresolved problem in X-ray astronomy for over two decades. We develop 
the model proposed Gicrlinski & Done (2004), which models the soft excess with 
heavily smeared, ionized, absorption, by including the emission that musthe associated 
with this absorption. We show that, rather than hindering the ionized absorption 
model, the addition of the emission actually helps this model reproduce the soft excess. 
The emission fills in some of the absorption trough, while preserving the sharp rise 
at ~1 keV, allowing the total model to reproduce the soft excess curvature from 
a considerably wider range of model parameters. We demonstrate that this model 
is capable of reproducing even the strongest soft X-ray excesses by fitting it to the 
XMM-Newton EPIC PN spectrum of PG1211-hl43, with good results. The addition of 
the emission reduces the column density required to fit these data by a factor ~2 and 
reduces the smearing velocity from ~0.28c to ~0.2c. Gierlihski & Done suggested a 
tentative origin for the absorption in the innermost, accelerating, region of an accretion 
disk wind, and we highlight the advantages of this interpretation in comparison to 
accretion disk reflection models of the soft excess. Associating this material with a 
wind off the accretion disk results in several separate problems however, namely, the 
radial nature, and the massive implied mass-loss rate, of the wind. We propose an 
origin in a 'failed wind', where the central X-ray source is strong enough to over- 
ionize the wind, removing the acceleration through line absorption before the material 
reaches escape velocity, allowing the material to fall back to the disk at larger radii. 

Key v^rords: galaxies: active - galaxies: Quasar - X-rays: galaxies - galaxies: 
PG1211+143. 



1 INTRODUCTION 



The origin of the soft X-ray (<2 keV) excess emission ob- 
served in many type-1 AGN has remained a long standing 
problem in X-ray astronomy despite considerable attention 
over the past two decades {e.g. Arnaud et al. 1985, Saxton 
et al. 1993, Mineshige et al. 2000). Several different ideas 
have been proposed for its origin including the thermal tail 
from the accretion disk {e.g. Pounds et al. 1986) and a sepa- 
rate cool Comptonizing region {e.g. Zdziarski et al. 1996) in 
either an optically thick transition region between the disk 
and a hot inner flow (Magdziarz et al. 1998) or in an ionized 
skin on the accretion disk (Janiuk, Czerny & Madejski 2001). 
These models encounter considerable problems; specifically, 
that the best-fit temperature for such thermal components 
is remarkably constant across a wide range of objects that 
exhibit a considerable range of black hole masses, luminosi- 



ties and accretion rates {e.g. Czerny et al. 2003; Gierlinski 
& Done 2004: hereafter GD04; Crummy et al. 2005). 

Detailed modelhng of the XMM-Newton EPIC PN spec- 
trum of the Narrow-line Seyfert 1 galaxy IH 0707-495 sug- 
gested that Compton refiection from a partially ionized skin 
on the accretion disk can reproduce the strong soft excess 
(Fabian et al. 2002). This model provides a natural explana- 
tion for the 'constant temperature' of the soft excess emis- 
sion since the 'excess' is now primarily associated with the 
abrupt change in opacity at ~0.7 keV. Partially ionized ma- 
terial is much more reflective below this energy, and the 
rising reflected continuum is enhanced by associated CNO 
lines and radiative recombination (RRC) emission features. 
Together, these can produce the soft excess, with the advan- 
tage that the characteristic energy of the emission is fixed 
by atomic physics, so is largely independent of the proper- 
ties of the AGN (Zycki & Czerny 1994). The narrow atomic 
features associated with ionized reflection are not evident in 
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the spectra of the soft excess suggesting that large (relativis- 
tic) velocity smearing is present in this material (Fabian at 
al. 2002). In addition, the largest soft excesses need a re- 
flected fraction considerably greater than unity, arguing for 
an unusual geometry such as i) a highly clumpy/corrugated 
disk surface so that the observed spectrum is dominated by 
multiple reflections of the mostly hidden intrinsic hard X-ray 
source (Fabian et al. 2004; Ross et al. 2002), or ii) a model in 
which strong light-bending towards the black hole increases 
the X-ray flux incident on the disk {e.g. Fabian ct al. 2005). 
The first of these models implies a complex physical con- 
flguration, perhaps as the result of disk instabilities that 
may be produced in high mass accretion rate sources such 
as Narrow-line Seyfcrt Is (NLSls), while the latter naturally 
produces the required large velocity smearing from the inner 
regions of the accretion disk. 

Gierliriski & Done (2004) proposed an alternative origin 
for the soft excess originating from similar atomic processes 
in partially ionized material, but with a different geometry. 
They suggested that the partially ionized material is seen in 
absorption, rather than in reflection. Again, a large veloc- 
ity dispersion is required to hide the atomic line and edge 
features, and results in a broad, smooth absorption trough 
that adds considerable curvature to the X-ray spectrum be- 
tween ~0.3-5 keV. The advantage of this model over reflec- 
tion is that a larger range soft excess sizes can be produced 
without invoking geometric changes. However, while the re- 
flection models include the self-consistent emission from the 
partially ionized material, the absorption model presented 
in GD04 does not. This is equivalent to assuming that the 
material subtends a small solid angle, in contrast with the 
large fraction of typc-1 AGN that show soft excess emission. 
Chevallier et al. (2005) include this emission as well as re- 
flection and suggest that all three components work together 
to produce the soft X-ray excess. 

GD04 suggest a tentative origin for the absorbing ma- 
terial in a UV line-driven accretion disk wind; that is mate- 
rial launched from the disk that is subsequently accelerated 
to high velocities by the radiation pressure on UV absorp- 
tion lines, in similar fashion to the winds from hot stars 
{e.g. Earners & Cassinelli 1999), reasoning that a large ve- 
locity dispersion is implicit to material launched from the 
inner disk and hence strong velocity smearing is to be ex- 
pected. The physical properties of such line-driven accretion 
disk winds have been investigated in considerable detail in 
other work {e.g. Arav, Li & Begelman 1994, Murray & Chi- 
ang 1995, 1997 & 1998, Murray et al. 1995, Proga, Stone & 
Kallmann 2000 - hereafter PSKOO, Proga & Kallman 2004 
- hereafter PK04, etc). The general properties implied for 
accretion disk winds match well to the blue-shifted (<10'^ 
kms~^) optical-UV- X-ray absorption lines seen in the Broad 
Absorption Line (BAL) quasars (Proga & Kallman 2004), 
and may also be responsible for the broad emission lines 
seen in all Quasars (Murray & Chiang 1997; Elvis 2000). 
Furthermore, such a wind provides an attractive and sim- 
ple way of dynamically linking the otherwise distinct regions 
within the standard AGN uniflcation scheme. 

In this paper we extend the absorption model proposed 
in GD04 to include the emission that must be associated 
with this absorption. In Section 2 we briefly present the 
details of the photoionization modelling and the velocity 
smearing. In Section 3 we examine the properties of the 



absorption and emission in more detail, focusing on the im- 
pact of the emission on the resulting spectral shape, and 
demonstrate that the combination of emission and absorp- 
tion removes many of the fine-tuning issues present in the 
pure absorption model of GD04. 

We test the impact of the addition of the emission on 
fits to real data by fitting the XMM-Newton spectrum of 
the bright, strong soft excess Quasar PG1211-I-143 with our 
models of the absorption and emission from a disk wind, 
together with ionized reflection and the narrow absorption 
systems known to be present in this source (Section 4). Sec- 
tions 5 and 6 discuss the physical origin of the absorbing 
material and its properties. 



2 MODELLING THE WIND 
2.1 Photoionized material 

We calculate the absorption and emission from photoionized 
material using the atomic code XSTAR* v2.1kn3 (Bautista 
& Kallman 2001). Given that we associate the smeared, 
blurred, ionized absorption with emission and absorption 
in the inner, accelerating, regions of a line-driven wind from 
the accretion disk, we chose gas conditions for the XSTAR 
model based in part on the values the data require and part 
on the simulations of AGN accretion disk winds presented 
in PSKOO & PK04. Specifically, the data require a column 
density range of Nh ^10^''*"^'' cm~^ (GD04) while the sim- 
ulations suggest a range of ionization states and a density 
of logio(C)=l-5-3.5 & n=10^^ cm^"^, respectively, for the in- 
ner regions of the wind. XSTAR self- consistently calculates 
the temperature of the material through the wind, we use 
solar metal abundances and we include the effects of mild 
turbulence in the gas (with a velocity dispersion of 100 km 
s~^). The parameters defining the ionizing continuum were 
chosen to best approximate the properties of the emission 
from the inner regions of a NLSl quasar. Specifically, the 
energy spectrum of ionizing X-rays is described by a power- 
law {i.e. N(E)=ifp;E~'") with a slope of r=2.4 (an energy 
spectral index of -1.4 in XSTAR) and an intrinsic luminosity 
of Li„t=10'"^ erg s"^ between 1-1000 Rydbergs. 

Absorption seen along a line-of-sight is (to first order) 
independent of the covering fraction of the absorbing mar 
terial, however the relative strength of the emission associ- 
ated with this absorbing material is sensitive to the covering 
fraction. We account for this fact later in this analysis and, 
initially, calculate both the absorption and emission spec- 
tra for a covering fraction corresponding to material cover- 
ing the entire sky from the point of view of the continuum 
source. A more detailed description of the treatment of cov- 
ering fraction is given in Section 3.3. 

There is considerable discussion in the community con- 
cerning the physical conditions under which photoionization 
calculations such as these should be calculated. Specifically, 
whether one can assume either constant density or constant 
pressure throughout the material. There is no clear consen- 
sus on which assumption best represents the real physical 



* http://heasaTc.gsfc.nasa.gov/docs/software/xstax/xstcir.html 
© 0000 RAS, MNRAS 000, 000-000 



Failed disk winds; a physical origin for the soft X-ray excess? 3 



conditions. Chevallier et al. (2005) argue for constant pres- 
sure conditions and highlight one of the more serious prob- 
lems with the constant density model presented in GD04, 
namely that this model is extremely sensitive to the free pa- 
rameters, particularly the ionization state, and that consid- 
erable fine tuning of these parameters is required in order to 
to match the observations of the soft X-ray excess. The con- 
stant pressure equilibrium models of Chevallier et al. avoid 
much of the fine tuning problem because the illuminating 
X-rays naturally induce a thermal instability within the 
material (e.g. Krolik, McKee & Tarter 1981). Beyond this, 
the material may fragment into small, cold, dense clumps 
with low covering factor (Rozahska et al. 2002, Rozariska et 
al. 2005, Chevallier et al. 2005). Thus the instability effec- 
tively leads to disruption of the lower ionization portion of 
the cloud, generating absorbed spectra with a similar overall 
shape and, crucially, a similar absorption trough depth for 
a wide range of ionization states and column density. 

Despite this advantage we choose to perform the pho- 
toionization calculations under the assumption of constant 
density primarily because the physical parameters are some- 
what simpler and more intuitive. In Section 3 we show that 
the addition of the emission component removes some of the 
required fine tuning in the constant density models, as the 
emission fills in part of the absorption trough. This effect 
becomes increasingly significant as the absorption trough 
gets deeper. This prevents a small change in the ionization 
state or the column density resulting in a large change in 
the absorption profile. 

The calculated spectra are dominated by OVII, OVIII 
and iron L-shell features seen in both emission (rest energies 
~0.55, 0.75 & 0.7-1.5 keV respectively) and absorption (rest 
energies ~0.74, 0.87 & 0.73-1.55 keV respectively). In ad- 
dition, the spectrum shows absorption associated with the 
iron M-shell unresolved transition array (~0.73-0.8 keV, Be- 
har, Sako & Kahn 2001). Figure 1, Upper Panel, shows an 
example emission and absorption spectrum from the grids 
of spectra calculated by XSTAR. 

2.2 Velocity structure 

Simulations of line-driven disk winds show a complex ve- 
locity structure (PK04), that is beyond the scope of this 
work to duplicate. The simulations demonstrate that the 
wind material typically accelerates in the radial direction, 
however there are several distinct velocity regions includ- 
ing a region in which the material in-falls, a low-density, 
rapidly accelerating region and a somewhat higher density, 
less-rapidly accelerating region. Further complicating mat- 
ters, while the final velocity-smeared absorption profile de- 
pends only on the velocity distribution along specific line 
of sight, the final emission smearing profile depends on the 
global velocity distribution of all the emitting material. For 
this work, we make the simplifying assumption that the dis- 
tribution of velocities present in the wind is a Gaussian, with 
a width that represents a typical velocity present in the wind 
(more precisely, 84% of the ionized material in the wind has 
a velocity given by the value, a^ind, or less). 

While this extremely simplistic assumption is unlikely 
to represent the correct profile in the case of absorption for a 
wind, as it results in as much red-shifted absorption as blue- 
shifted absorption, it may be more reasonable in the case 
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Figure 1. Upper Panel: An example of the emission (blue) and 
absorption (red) X-ray spectra from gas photoionized by the X- 
ray continuum of a quasar. The intrinsic input spectrum (black, 
dotted) is a power-law with V=2A & Kp;=10~^ photons keV~^ 
cm~^ s~^. The photoionized material has Nh=10^'' cm~^ & 
loglo(?)=2.6, in addition to the details given in Section 2.1. Lower 
Panel: The photoionized emission (blue) and absorption (red) 
spectra from the upper panel, smeared with a Gaussian velocity 
distribution with a typical v/c=0.2 The intrinsic input spectrum 
(black, dotted) is the same as in the upper Panel. 

of the total emission profile which will necessarily include 
emission from all parts of the wind. We note that changing 
the global velocity profile will result in significant differences 
in the absorption line profiles and, to a lesser extent, the 
emission line profiles. 

Figure 1, Lower Panel, shows the effect of convolving 
the velocity distribution with the example XSTAR absorp- 
tion and emission spectra and demonstrates that the in- 
trinsically narrow features appear as continuum-like com- 
ponents if the velocity distribution is sufficiently broad. In- 
terestingly, between 0.7 and 1.3 keV, the smeared emission 
spectrum rises somewhat more steeply than the curvature 
imprinted purely by absorption, suggesting that the emis- 
sion will strongly affect the shape of the soft X-ray excess 
predicted by this model, provided the covering fraction is 
sufficiently large. 

We note that this model is not the same as one in which 
the material has a large turbulent velocity. The amount of 
turbulent velocity present in absorbing material strongly af- 
fects the opacity of heavily saturated absorption lines (see 
e.g. Peterson 1997). A large turbulent velocity will broaden 
the saturated cores of these lines resulting in considerably 
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more absorption than we might expect from this material. 
The impact of turbulence is discussed in more detail in Sec- 
tion 5. 



3 A DETAILED EXAMINATION OF THE 
WIND PROPERTIES 

Before applying the disk wind model it is instructive to un- 
derstand how the model varies with the free parameters, in 
particular, how the column density (N^f), the ionization pa- 
rameter (^) and the global covering fraction of the material 
(C/ - see Section 3.3) affect the emission spectrum and how 
significant this component is to the resulting total spectrum. 
The total emitted luminosity is given by; 



= Cf {LrrC + Llines + Lff) 



(1) 



where Lums is the luminosity of the emission (bound-bound 
transitions), Lrrc is the luminosity of the radiative recom- 
bination continua (free-bound transitions) and L// is the 
luminosity of the free-free continuum (free-free transitions) . 
All three emission mechanisms have the same dependence on 
the electron density of the emitting material (rie) and the 
volume of the emitting material {V), which we can relate 
directly to the free parameters in out wind model. Specifi- 
cally, 

Lrrc = nlV FRRc{^,hv,aRRc) 

Llines ^ ^lines (d Cannes') 

Lff = nlVFff{^,hu,aff) (2) 

where Frrc, Funes &i Fff are functions of the ionization 
state of the material, f , the energy released by the interac- 
tions involved, hi/, and a coefficient appropriate to each type 
of interaction, a,,, (see Appendix A). 

For the XSTAR models used here, the volume of the 
emitting material is V—A-jir^ AR, where r is the radius from 
the central source to the inner edge of the wind material and 
AR is the radial depth of the shell. In addition, A_R can be 
defined in terms of the density and the column density of 
the material such that AR — Nn/n. Combining these with 
Equations 1 & 2 we have. 



Cf rie ^irr Nh (Frrc + Fu 



F, 



ff) 



(3) 



Furthermore, the ionization state of the material is 
given by ^ = Li„t /ner^ , where Li„t is the intrinsic lumi- 
nosity of the ionizing continuum. Combining this with Eq 3 
gives; 

_ CfN„ 



Lint 



■4,-rv{FRRc: + Fii„es + Fff) 



(4) 



showing that, to zeroth order, the luminosity in the emis- 
sion, as a fraction of the intrinsic ionizing luminosity, is pro- 
portional to the column density and the covering fraction 
and inversely proportional to the ionization state. We em- 
phasise that the dependence on ionization state is more com- 
plex than the simple inverse relation implied by the initial 
part of Eq 4, in that the possible atomic transitions are dif- 
ferent for different ionization states of the gas, giving very 
different values for each of the Fn-functions as we vary 5. 
Building on the simple relationships in Eq 4, we examine 
the impact of, and the interplay between, the emission and 




Energy (keV) 

Figure 2. The effects of varying N^f on the smeared wind model. 
Upper Panel: The dotted line represents the intrinsic input con- 
tinuum (power-law, r=2.4, Kp;=10~^ photons keV~^ cm~^ 
The solid lines show the smeared absorption and emission spec- 



tra (logio (0=2.6, o-„i„d=0.2, 
3.3 & 6.6 Xl022 (black, red 



Cf=1.0) for column densities of 1, 
: green), 1, 3.3 & 6.6 xlO^^ (dark 
blue, light blue & purple) & 1x10^"* (brown) cm~^. Lower Panel: 
The solid lines show the 'total' spectrum (absorption X intrinsic 
spectrum -|- emission) for the same set of parameters. 



absorption on the shape of the total X-ray spectrum, as we 
vary the N_fr, ^ and Cf. 



3.1 Column density 

Figure 2 shows the spectral shape changes associated with 
changing the column density. The upper panel shows the ef- 
fect of changing Ni^ on the absorption and emission spectra, 
while the lower panel shows the effect on the total spectrum. 
Below Nh~5x10^^ cm~^, even with a covering fraction of 
one, the smeared emission is not a significant correction to 
the shape of the absorbed spectrum. Above this however, the 
emission becomes increasingly important, eventually becom- 
ing the dominant spectral component between 0.2-3 keV at 
column densities of ^10^"* cm~^. 

The absorption increases dramatically with increasing 
column, as expected, but there is also a more subtle ef- 
fect, namely that the trough of the absorption component 
shifts to lower energies with increasing column, indicating 
the presence of lower ionization material. This is a natural 
consequence of the increase in absorption in the front layers 
of the cloud reducing the flux transmitted deeper within the 
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cloud (the reduction in flux due to the increasing distance 
is neghgible since the cloud has AR<<R). Thus increasing 
the column leads to more shielding and so to progressively 
lower ionization material at the back of the cloud. 

The emission also follows this pattern, with the emit- 
ted flux increasing strongly with column, but also showing a 
subtle shift to slightly lower ionization states, as revealed by 
the slightly lower peak energy (due to the strong OVII/OVIII 
lines) Interestingly, the emission is always characterised by 
a sharp rise below ~1 keV, precisely the point where we ob- 
serve the sharp rise of the soft X-ray excess. This suggests 
that the emission component may, in fact, assist us in ac- 
counting for the soft X-ray excess! Intruigingly, similar mod- 
els of emission (though not absorption) from photoionised 
material were used to fit the soft excess seen in the X-ray 
spectra of many accreting pulsars (Hickox, Narayan & Kail- 
man 2004). 

The co-added emission plus absorption spectra shown 
in the lower panel of Figure 2 support this, showing that 
the infill from the emission offsets the effect of lower ioniza- 
tion material at high columns, leading to an apparent rise 
in the spectrum which is always close to ~1 keV. Signifi- 
cantly, however, the peak in the emission is quite obvious 
in these model spectra, leading to a pronounced dip below 
~0.6 keV. This should be an observable consequence of a 
large contribution from the emission of the material. While 
such a dip is not commonly reported in the shape of the soft 
X-ray excess in AGN, this could be masked by the presence 
of either a Compton reflection continuum, narrow ionized 
absorption intrinsic to the AGN, neutral absorption intrin- 
sic to the AGN and/or galactic absorption (typically IQ^""^^ 
cm-2). 

3.2 Ionization state 

Figure 3 shows the spectral shape changes associated with 
changing the ionization state. The upper panel shows the 
effect of changing ^ on the absorption and emission spectra, 
while the lower panel shows the effect on the total spectrum. 
The effects of increasing ^ are broadly opposite to the effects 
of decreasing the column density (as we would expect from 
equation 4), in that the absorption trough becomes less pro- 
nounced and shifts to higher energies. Physically, the opac- 
ity of ionized material falls with increasing ionization state 
and the material becomes dominated by higher ionization 
species of each ion, giving rise to less pronounced, higher 
energy, absorption troughs at higher ionization states. 

While increasing the column led to some lower ioniza- 
tion material present at the back of the cloud, this is a very 
small effect compared to lowering the ionization parameter 
at the front of the cloud. Fig. 3 shows that the low energy 
absorption increases substantially below logio^ ~2.1, mainly 
due to the presence of L shell opacity from Oxygen VI and 
Nitrogen V as well as Fe M shell features. 

The effect of varying ^ on the emission is considerably 
more 'messy' than varying Nh- Specifically, the contribution 
of N VII, OVII & OVIII, relative to highly ionized species of 
sodium, magnesium and silicon, drops dramatically above 
ionization states of logio (5)~2.8. This removes the peak in 
the emission at ~0.6 keV due to OVII line emission, so the 
spectrum is much flatter at these ionization states. 

The total spectrum (seen in the lower panel of Fig. 3) 
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Figure 3. The effects of varying ^ on the smeared wind model. Up- 
per Panel: The dotted line represents the same intrinsic input con- 
tinuum detailed in Fig 2. The solid lines show the smeared absorp- 
tion and emission spectra (N^=10'^^ cm~^, cr^i„ii=0.2, Cj=1.0) 
with ionization states of Iogio5=1.5, 1.83, 2.16, 2.5, 2.83, 3.16 
& 3.5 (black, red, green, dark blue, light blue, purple & brown). 
Lower Panel: The solid lines show the total spectrum for the same 
set of parameters. 

shows again that including the emission counteracts some 
of the sensitivity to the parameters predicted from absorp- 
tion alone. There is a soft excess which rises at '^l keV for 
logio (5) = 1.8-3. Again, a key characteristic of the emission 
from ionization states below logio(5)~2.6 is that it peaks at 
~0.6 keV, and these models predict a dip at lower energies 
if the emission is important. 

Including the self-consistent emission from the partially 
ionized material results in a stable soft excess for a wider 
range of logio(^) than expected from pure absorption alone. 
We note that the slope of the smeared emission spectrum 
between 0.7-2 keV is relatively sensitive to changes in ioniza- 
tion state which may explain the distribution of soft excess 
shapes observed in AGN (Crummy et al. 2005). 

3.3 Covering fraction 

As discussed in Section 2.1, the properties of the absorption 
seen along the line-of-sight (l.o.s) are, to flrst order, inde- 
pendent of the covering fraction of the material, however 
the relative strength of the emission is sensitive to the cov- 
ering fraction. In fact, the absorption spectrum does depend 
weakly on the local covering fraction (i.e., whether the ma- 
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terial in the immediate environment of the Lo.s is clumpy 
or uniformly smooth) due to some subtleties in the way XS- 
TAR propagates radiation from one layer in the material to 
the next (Kallman priv. comm.). By contrast, the emission 
is sensitive to both the global covering fraction (i.e., whether 
the global distribution of the material is in a uniform shell, 
a non-uniform shell, a conical geometry, etc) and the local 
covering fraction. The covering fraction variable available in 
XSTAR (cfrac) strictly refers to the local covering fraction. 
The global covering fraction is, by default, unity for a sin- 
gle XSTAR run, since XSTAR simulates a spherical shell of 
material, however the grid of calculated emission and ab- 
sorption spectra are stored as separate table model and are 
treated individually by XSPEC. In this case, the global cov- 
ering fraction is folded into the normalization of the emission 
additive model. 



3. 3. 1 The normalization of the emission 

It is instructive to understand what the normalization of 
the emission component actually means before we use it in 
XSPEC. The normalization of the emission {Kemm) is given 
by; 



Cf Lint 1 



1038 



kpc 



(5) 



where C/ is the global covering fraction and Dtpc is the 
distance to the source in kpc (this expression is derived in 
detail in the XSTAR 2.1kn3 Manual, § 7.2.3.1"). For a power- 
law continuum of the form N{E)=Kpi{E /IkeV)'^ , where E 
is energy in keV and Kpi is the power-law normalization at 
IkeV in photons keV~^ cm~^ s~^, Li„t is given by; 
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Figure 4. The effects of varying C f on the smeared wind model. 
Upper Panel: The dotted line represents the same intrinsic input 
continuum detailed in Fig 2. The solid lines show the smeared 
absorption and emission spectra (Nj:/=10^^ cm~^, logio (§)=2.6, 
(T„=0.2) with covering fractions of 1.0, 0.85, 0.7, 0.55, 0.4. 0.25 
& 0.1 (black, red, green, dark blue, light blue, purple & brown). 
Lower Panel: The solid lines show the total spectrum for the same 
set of parameters. 



U^t = Dip, Kpi ^ ^ (_r + 2) 

where the constant X contains the conversion from cm^ to 
kpc^ and from keV to erg (=1.53x 10^"), E2 & E\ are the up- 
per and lower limits of the energy range respectively (1-1000 
Rydbergs in XSTAR). Combining Eq 5 & 6, the distance de- 
pendence of the emission normalization can be removed and 
we obtain an equation that depends solely on the details of 
the input power-law continuum (specifically its normaliza- 
tion and its slope) and the global covering fraction of the 
ionized material. 

Clearly, allowing the normalization of the emission to 
be independent of the normalization of the underlying con- 
tinuum in any actual fit to real data, gives us a measure of 
the global covering fraction of the ionized material. In the 
remainder of this work, we assume that the ionized material 
is uniformly smooth on all scales and interpret the global 
covering fraction as an indicator of the large scale geometry 
of the wind material. 



3.3.2 The effect of changing Cf 

Figure 4 shows the spectral shape changes associated with 
changing the covering fraction of the ionized material. The 
upper panel shows the effect of changing C/ on the emission 
component (the absorption component being independent 
of the global covering fraction) , while the lower panel shows 
the effect on the total spectrum. The emission has a consid- 
erable effect on the spectrum even at relatively low covering 
fractions, particularly on the depth of the total observed ab- 
sorption trough, however it is only at covering fractions of 
>50% that the emission begins to dominate the spectrum. 
The characteristic 'double-humped' absorption profile be- 
comes considerably less pronounced as the emission begins 
to fill in the absorption at low energies. We note, again, that 
regardless of the relative strength of the emission the spec- 
trum in the 0.7-1.5 keV range shows a steep rise reminiscent 
of the soft excess. Again, when the emission becomes very 
strong the spectrum shows an increasingly pronounced dip 
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in the total spectrum below ~0.6 keV, similar to the dips 
discussed previously for Nh & ^• 



4 PG1211+143; A CASE STUDY 

The objects which will provide us with the tightest con- 
straints on the wind parameters are the bright quasars with 
particularly large column densities. On this basis we choose 
PG1211-I-143 on which to tost this model. 

PG1211-I-143 is a bright (V,„a9=14.38), nearby 
(z=0.089) quasar that has been identified as a Narrow Line 
Seyfert 1. The X-ray spectrum of PG1211-I-143 is known to 
be complex. RXTE observations suggest a significant Comp- 
ton reflection component (Janiuk et al. 2001). ASCA obser- 
vations revealed a strong soft X-ray excess and evidence that 
the iron Ka emission line was broad (Reeves ct al. 1997). 
XMM-Newton observations revealed the presence of several 
complex absorption systems that can either be interpreted 
as high velocity (-^24,000 km s~^), highly ionized absorption 
systems (Pounds et al. 2003), or lower velocity (~3, 000-6, 000 
km s~^) systems with a somewhat lower ionization state 
(Kaspi et al. 2005). 

Most importantly, PG1211-I-143 has one of the strongest 
soft X-ray excesses observed, with the X-ray flux below 0.6 
keV being up four times greater than the flux predicted 
by continuum fits to the 1-10 keV spectrum (Pounds et 
al. 2003). Such a large soft excess makes it the ideal test 
of models for the origin of this component and as such it 
was used by GD04 as a test of their absorption-only model. 
Here we extend their fits to include the smeared emission as 
well as absorption, testing the ability of the full model to fit 
the soft excess seen in this object. 

4.1 Data Reduction 

XMM-Newton observed PG1211-I-143 in 2001 for 55ks (Obs 
ID: 0112610101). During the observation the EPIC PN was 
operated in Large Window Mode with the medium filter. 
The EPIC PN data were reduced and processed with the 
Science Analysis System (SAS v6.5.0) using the standard 
processing chain and the most up to data calibration data 
available as of November 2005. Periods of high background 
were excluded from the analysis and the spectrum was ex- 
tracted from a 45" extraction region resulting in a net expo- 
sure time of ~49.5 ks. The spectrum was background sub- 
tracted and binned to have a minimum of 20 counts per bin 
to allow the use of chi-squared fitting statistics. 

4.2 Fitting the X-ray spectrum 

We begin by defining a model which accounts for all the 
complexity known to be present in the 0.3-12.0 keV spec- 
trum of PG1211-I-143. We use a power-law continuum mod- 
ified by Compton reflection, two un-smeared, partially ion- 
ized absorption systems, together with Galactic absorption 
(2.85x10^° cm-2 - Murphy et al. 1996). 

It is unclear from either Pounds et al. (2003) or Kaspi 
ct al. (2006) which of their interpretations for the nature 
absorption complexities best represents the data. However, 
given that our aim is to test whether the accretion disk wind 
model can represent the broad shape of the X-ray spectrum. 



the details of the complex absorption systems are not par- 
ticularly relevant to the work presented here, provided that 
they are sufficiently well modelled that they do not damage 
the quality of the overall fit. To this end, we adopt the inter- 
pretation of Pounds et al. on the basis of simplicity (since 
this method does not require ion by ion fitting). We con- 
struct two sets of XSTAR table models covering the param- 
eter ranges detailed in Pounds et al. (Nif,i=5xl0^^ cm~^, 
^1 ~3.4, Nh,2=6x10^^ cm"^, ^2 ~1.7). We include the emis- 
sion associated with each of those narrow absorption systems 
and allow the normalization of each of the emission compo- 
nents (equivalent to the covering factor of the material) to 
be a free parameter in the model. 

The impact of Compton reflection is strongly dependant 
on the viewing angle. In the standard unified AGN picture, 
the accretion disk is seen almost face on in type-1 AGN 
such as PG1211-I-143, while the the region of the (putative) 
torus exposed to the X-ray continuum emission is almost 
edge-on to our line of sight. Thus we expect the total re- 
fiection signature to be dominated by the component asso- 
ciated with the accretion disk, which may be ionized (Januk 
et al. 2001), and should be relativistically smeared by the 
large orbital velocities. We use the XION ionized reflection 
model (Nayakshin et al. 2001) to represent this component. 
This is an extremely flexible and sophisticated model, incor- 
porating several different geometries for the accretion disk 
and the X-ray source(s), variable iron abundance and a wide 
range of accretion rates. We choose a set of fixed values for 
the majority of the model parameters that are appropriate 
for the properties of this source. Specifically, we choose a 
magnetic flare based geometry for the X-ray source (which 
forces f2/27r=l), we fix the outer edge of the accretion disk 
at 1000 Rs , the height of the flares above the accretion disk 
at 5 Rs, the accretion rate to be Q.l&MEdd (based on an es- 
timate oi'Lboi/^Edd)^ , the Fe abundance to be solar and the 
viewing angle of the disk to be almost face-on (cos(i)=0.85). 
The remaining free parameters in the fit are; the radius of 
the innermost edge of the accretion disk (this governs the 
amount of relativistic broadening observed) and the ratio of 
the X-ray source luminosity to the luminosity of the disk 
(Li^/Ldisfc). 

On top of this base model we add the absorption and, 
crucially, the emission from the disk wind model discussed in 
Section 2. The ionization state, column density, the velocity 
smearing width and the normalization of the emission are 
free parameters in the model. 

Figure 5 shows the best-fit model continuum compo- 
nents (top panel), the total best-fit model (middle panel), 
the model fit to the EPIC PN data (bottom panel, up- 
per section) and the data/model ratio residuals (bottom 
panel, lower section). Table 1 details the best-fit parame- 
ters of the model. The model provides a good fit to the 
XMM-Newton EPIC PN data (xLd=l-08). The curvature 
introduced into the spectrum from the smeared absorp- 
tion/emission model does an excellent job of reproducing 
the apparent soft X-ray excess. The only remaining residu- 
als are narrow features around 2.2 keV (a known calibration 



t Ljo;=4.4xlO*s erg s^^ & LBdd=5.7xlO*s erg s"! for 
MBif,pGi2ii=4.5xlO^ M0; Kaspi et al. 2000 
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Figure 5. Top Panel The best-fit model continuum components, 
namely the underlying power-law continuum (dashed black line) 
modified by the smeared absorption (red), refiection from the ion- 
ized skin of the accretion disk (green) and the emission associated 
with the smeared wind system (blue). Middle Panel The complete 
best-fit model fit to PG121H-143 including both narrow ionized 
absorption systems, and the emission associated with this absorp- 
tion, and galactic absorption. Bottom Panel The fit of this model 
(solid line) to the EPIC FN spectrum of PG1211+143 (upper sec- 
tion) and the data/model ratio residuals (lower section). 

issue around the gold edge in the EPIC PN response^, and 
~7 keV, which is probably due to the very highly ionized ab- 
sorption system (Pounds et al. 2003) not being completely 

§ Document 

ref: XMM-SOC-CAL-TN-0061. http://xmm.vilspa.esa.es/ exter- 
nal / xmm^w_cal /calib/ documentation / 



described by our model. We speculate that this is due to the 
incomplete Fe K-shell physics in the version of XSTAR used 
to perform the photoionization calculations here (v2.1kn3 
as opposed to the newest version available, v2.11, which in- 
cludes more detail for Fe K-shell transitions: Kallman et 
al. 2004). 

We note that our spectra decomposition attributes the 
majority of the soft excess to smeared absorption and emis- 
sion rather than with the reflection component. Signatures 
of the reflection are clearly present in the spectrum (see 
Section 4.4), however the best-fit parameters of the reflec- 
tion continuum are such that this component does not rise 
strongly below 0.7 keV. A note of caution is required here 
however because both refiection dominated and absorption 
dominated models are both capable of providing a good fit 
to these data. This implies that either the models (in par- 
ticular the details of the ionization and velocity structure of 
the wind and/or the accretion disk) or the data are of in- 
sufficient quality to distinguish these two physics scenarios 
unambiguously (Sobolewska & Done 2006). 

4.3 The smeared absorption/emission system 

Fig. 5 shows that the apparent soft excess is well modelled 
by the combination of the smeared absorption and emis- 
sion from a column with N£f~1.5xl0^^ cm~^ and velocity 
smearing of ~0.2c. The solid angle implied by the amount 
of emission is ~0.4, so a substantial fraction of the central 
source is covered by this material, as required by the ob- 
served ubiquity of the soft excess. The pronounced dip be- 
low ^O.Q keV predicted for a smeared absorption/emission 
system with this covering fraction is not obvious in the data, 
however the model provides a good fit, demonstrating that 
this feature is masked in the observed spectrum by the re- 
fiection component and the impact of galactic absorption. 

Comparing these best-fit parameters with those for the 
absorption- only model presented in GD04, highlights the im- 
portant effect that the emission has on the implied physical 
properties of the wind. The required column is reduced by 
a factor 2 by including the emission (from 3.3 x 10^^ cm~'^ 
in GD04), though its ionization state is not significantly af- 
fected (560 here c.f. 460 erg cm s"^ in GD04). More im- 
portantly, the velocity smearing is reduced by a factor 30% 
(from 0.28c in GD04). This is primarily the result of the 
emission component helping to fill-in the strong absorption 
around 0.9 keV, reducing the level of smearing required to 
blend the two strongest absorption regions together. Thus 
including the emission makes the wind significantly less ex- 
treme in column and velocity than inferred from using the 
absorption alone. 

4.4 Compton refiection 

In line with what we expect from Narrow-line Seyfert 1 
objects, the underlying continuum emission is steep, with 
r~2.4, and the spectrum displays clear signatures at- 
tributable to Compton refiection from the accretion disk. 
The iron Via line is broader than the instrument response, 
but not to the level expected from emission from the last sta- 
ble orbit even from a Schwarzschild black hole. The lower 
bound on the disk radius for the line is 2QRg , so does not re- 
quire extreme spin, and the model itself by construction has 
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Table 1. Best-fit spectral parameters 
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0/27r = 1. Both these results are, however, model depen- 
dent. These data can be equally well fit without smeared 
absorption/emission using a continuum which is reflection 
dominated and requires extreme spin (Sobolcwska & Done 
2006). Plainly the origin of the soft excess has profound im- 
plications for decoding the geometry of the space-time in 
the vicinity of a super-massive black hole, as well as for the 
details of hard X-ray source(s) and the accretion disk. 

4.5 The narrow absorption systems 

The origin of the soft excess also has implications for 
the narrow absorption systems. Both the smeared absorp- 
tion/emission model and the smeared ionized reflection 
model contain an intrinsically stoop continuum spectrum 
(typically modelled with a power-law with r'~2. 2-2.4), while 
interpreting the soft excess as a separate continuum com- 
ponent (as the tail of the multi-colour blackbody emission 
from the accretion disk, for example) suggests a much flatter 
spectrum, (r~1.8-1.9). The ion populations of the narrow 
absorption systems depend strongly on the spectral index of 
the ionizing continuum and will yield subtly different narrow 
absorption systems in each case. Both Pounds et al. (2003) 
and Kaspi et al. (2006) interpret the soft excess as a sepa- 
rate continuum component and, as a result, the properties 
of the narrow absorption systems fit here are different to 
those that Pounds et al. & Kaspi et al. derive from the same 
data. 

In particular, while the ionization states and column 
densities found here arc rather similar to those of Pounds et 
al. (log(0-^3.27 & 1.99 c.f. ~3.4 & 1.7 and N^^lxlO^^ & 
1x10^2 c.f. ~5xl0^^ & 0.6x10^^ cm-2), the inferred outflow 
velocities arc different (45,000 & 13,000 km/s c.f. 24,000 & 
24,000 km/s, respectively). The steeper continuum in our 
model prevents a strong FeXXVI Lya (6.96 keV) transition 
being present in the spectrum at logig ~ 3.3. Instead, the 
spectral feature at ~7 keV (observed frame) is modelled with 
the wealth of Fc absorption lines present in our model at a 
rest frame energy of ~6.6 keV, resulting in a considerably 
higher apparent outflow velocity for the most ionized narrow 
absorption system. Again, we note that the outflow velocity 



could change using the newest K shell calculations (Kallman 
et al. 2004). 

None of the narrow absorption systems has significant 
emission detected from it. At CCD resolution the individual 
lines cannot be resolved, but we can place fairly stringent 
upper limits on the covering fraction of this material at less 
than ~6% of the sky. Examining the emission spectra in 
detail reveals that the covering fraction of the more strongly 
ionized gas is constrained primarilly by the lack of strong 
ionized iron K features around ~7 kcV, whereas the covering 
fraction of the more weakly ionized material is constrained 
by the lack of strong soft X-ray features, particularly those 
related to OVII (~0.7 keV) and OVIII (~0.9 kcV). 

We note that the turbulent velocity used in the pho- 
toionization modelling of the narrow absorption systems 
does aifect the observability of the emission spectrum, how- 
ever the effect is weak with velocities up to ~10,000 kms~^ 
resulting in a factor <2 change in the upper limit. 



5 CONCERNING DISK WIND PROPERTIES 

Given our best-fit wind properties we can make a crude es- 
timate of the inferred mass-loss rate, Mout. Assuming that 

our velocity smearing represents a geometrically thick (i.e. 
Rout»Rin) outflowing wind, the mass-loss rate approxi- 
mates to; 

Mout = Cf ATvRin Uenip V (8) 

where C/ is the global covering fraction of the wind, Rin 
is the wind launch radius, rie is the electron density, rup 
is the mass of the proton and v is the average velocity 
of the wind material. Using the parameters for our pho- 
toionized material (Rin=5xl0^'' cm and n ^lO^'^ cm^"^), 
Afout=190 M0 yr~^. By contrast, the mass accretion rate 
which corresponds to accretion at the Eddington limit for 
PG1211-I-143 (assuming an accretion efficiency of ~6%) is 
~1.3 Mq yr~ , presenting us with the uncomfortable situa- 
tion that the disk is losing over 100 times more mass through 
the wind than is being accreted. 

Another way to see the same conflict is via the column 
density, N_ff , through the wind, 

/•Rout pRout jTr 

N„= n,dR= dR (9) 

For a constant velocity outflow extending to cxd then 
N/f=1.2xl0^^ cm~^, dramatically larger than the measured 
column needed to produce the soft excess (~10^^ cm~^). 

Clearly the main problems lie in the assumption of a 
uniform density of the material and that the wind extends 
to 00 with a single, large velocity. Clearly a more detailed 
description of the density structure is required in order to 
evaluate the mass-loss rate from such a wind more accurately 
{e.g. the description given by Murray & Chiang - 1995, see 
Section 5.1), however if the wind slows down, as we might 
expect given that it has to escape from a large gravitational 
potential well, then the predicted M and Nh increase even 
further! 

Chevallier et al. (2005) perform a similar set of calcu- 
lations for their pure- absorption, constant pressure, models 
(using a standard set of parameter values) and again they 
calculate a mass-loss rate of Mout^W^ Msdd unless they 
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Table 2. Comparing the physical parameters of PG1211+143 
with those of accretion disk wind simulations (PK04) 
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incorporate an extremely low volume-filling factor. Consid- 
ering that we derive our wind parameters from detailed spec- 
tral fitting, and our model includes emission as well as ab- 
sorption, we consider the agreement between our results and 
the Chcvallior et al. calculations to bo rather good. This 
presents a major problem for the validity of interpreting 
and origin for the smeared absorption/emission as arising in 
a disk wind. 

Examining the model in detail a possible method of cir- 
cumventing much of this problem suggests itself. As noted in 
Section 2.2, the amount of turbulent velocity present in the 
material strongly affects the opacity of heavily saturated ab- 
sorption lines and analysing such a spectrum with a model 
that only incorporates a small turbulent velocity will result 
in a considerable overestimate of the column density of the 
material. Simulated XSTAR spectra with turbulent veloci- 
ties of ~3xl0* km s~^ (the maximum allowed in the code) 
show that this effect can result in up to ^3-4 times more 
absorption in the 0.3-3 keV band, for gas with similar prop- 
erties to the gas we use in our simulations, than for gas with 
small turbulent velocities. 

Fortuitously, whilst the velocity dispersion an accre- 
tion disk wind will differ considerably from that imposed by 
turbulent velocity, the extremely simple Gaussian smearing 
used in this work is actually rather similar to the velocity dis- 
persion from turbulent velocity (although, obviously, it does 
not incorporate the increase in opacity for heavily saturated 
lines). The required velocity shear we apply (^6x10'' km 
s^^), if interpreted as a turbulent velocity instead of a wind 
outflow velocity, might well be causing us to overestimate 
the column density of our material, and hence out mass-loss 
rate, by up to a factor ~10. Clearly this correction is not car 
pable of producing the three orders of magnitude required 
to solve the mass-loss rate problem. 

5.1 Comparison with simulations 

Detailed hydrodynamic simulations are a relatively new ad- 
dition to the study of accretion disk winds, however a com- 
parison between the spatial simulations presented in PSKOO 
& PK04, and spectral simulations presented here, is instruc- 
tive. The basic properties of the spatial simulations are rel- 
atively well matched to PG1211-I-143, with only the ratio 
^x/l^disk and the spectral index of the ionizing X-ray radi- 
ation being significantly different (sec Table 2). 

The spatial disk wind simulations show a wind with an 
extremely complex density and velocity structure, clearly 
showing that the constant density, thin shell, treatment as- 
sumed in our photoionization calculations is far too sim- 
plistic. Not only is the density structure complex, but the 
typical density of the the simulated wind is rrmch lower than 
the value assumed during our modelling (lO'' cm'^ c.f. 10^^ 



cm^"^), except in the very innermost regions. We note, how- 
ever, that the harder spectral index of the spatial simulations 
will result in a more strongly ionized wind, driving down the 
density of the wind in comparison with the values expected 
from a considerably softer X-ray source. The spatial simu- 
lations also show that the global wind velocity structure is 
extremely complex and strongly time dependant, in strong 
contrast with the simplistic velocity treatments discussed in 
Section 2.2. 

The region of the spatial simulations that best matches 
our ionization and column density conditions (logio(0~2-4 
and/or a column of Nh^W^^ cm~^) is the fast, relatively 
geometrically thin, high-density 'stream' which, in addition 
to having approximately the correct ionization state and col- 
umn density, also displays the required large velocity disper- 
sion and peaJis at speeds of ~30,000 km s~^ (this is, in fact, 
a lower limit on the wind velocity since there is a numerical 
ceiling on the velocities allowed in the sirimlations) . Wc note 
that despite this apparent match, for the particular simula- 
tion shown in PK04, there is no specific l.o.s through the fast 
stream which has both the observed column and the correct 
ionization state simultaneously. While the global wind veloc- 
ity structure in the spatial sirrmlations is extremely complex 
and strongly time dependant, making the correct velocity 
smearing for the emission difficult to calculate, absorption 
in the material in the fast stream is a simpler case. 

Murray et al. (1995) and Murray & Chiang (1995) de- 
scribe the velocity and density structure of an accretion disk 
wind along a single streamline. The fast stream in the spatial 
simulations is extremely similar to this case, allowing us to 
combine the Murray & Chiang descriptions of the density 
and velocity distribution along a streamline into a single 
density-weighted radial velocity smearing function (which 
turns out to be proportional to R~^) for the fast stream. 
Applying this smearing function in future work will be a 
considerable improvement in the treatment of absorption by 
accretion disk winds in AGN and, provided that the emission 
from this fast stream dominants the emission from the wind 
material as a whole, may also be appropriate for the ab- 
sorption. A more fundamental problem here is that many of 
the PG Quasars show soft X-ray excesses that arc very sim- 
ilar to that of PG1211-I-143. In the interpretation as a disk 
wind, this implies that all these sources see through a similar 
region of the wind, exhibiting a similar amount of absorp- 
tion in material with a similar ionization state. While the 
covering fraction of the fast stream is of order ~0.2 (PSKOO, 
PK04), the section which closely matches our ionization and 
column density conditions is much smaller. Unfortunately, 
while the fast stream of a line-driven disk wind may have 
some of the properties required to reproduce the soft X- 
ray excess in PG1211-I-143, it is unlikely that the l.o.s to so 
many PG Quasars would fall within the tightly constrained 
range of solid angle covered by the fast stream of such a 
wind. The situation becomes even more problematic when 
we consider that all these objects are Type-1 AGN, implying 
that the l.o.s is close to the normal to the plane of the disk. 
In contrast, the fast stream in the disk wind simulations is 
generally quite equatorial, rising no more than 25° from the 
disk plane. The simulations of line-driven winds seem to bear 
a striking resemblance to the outflows required in broad ab- 
sorption line quasars (PSKOO, PK04), but not to the rather 
faster, denser and larger scale height material required here. 
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In addition, the total mass-loss rate in the spatial 
simulations is only 0.5 M© yr~^, in stark contrast to the 
190 Mq yr~^ predicted by out spectral model. At first 
glance, it seems that the huge mass outflow rates implied 
by our spectral modelling can be simply reduced by drop- 
ping the density used in the photoionization calculations to 
a density more typical of that seen in the wind simulations. 
Reducing the density used in the spectral simulations to 
'^lO^ cm~^ would result in a reduction of the mass loss rate 
to Moutfiow^O.1-1 Mq yr^^. We then require, however, that 
the material is distributed in a thick shell, with AR>>R in 
order to obtain the required column density of ~10^^ cm~^. 
Unfortunately, relaxing the 'thin shell' condition generates 
its own obstacles, namely i) a considerable increase the com- 
putational time required to calculate the grid models, and ii) 
photon- loss, and hence energy-loss, through the upper/lower 
'surface' of the wind (see Appendix B & King & Pounds 
12003 for a brief description of this problem). 

Once again, however, we note that the more strongly 
ionized wind that results from the harder spectral index of 
the spatial simulations, not only drives down the density 
of the wind, but also leads to less material reaching escape 
velocity before becoming over-ionized. As the material be- 
comes increasingly more ionized, it becomes transparent to 
the UV line emission (as there are few ion stages left with 
substantial line opacity), reducing the mass- loss rate from 
such a wind in comparison with a wind exposed to a softer 
X-ray continuum. 

In fact, a sufficiently powerful hard X-ray flux causes 
the wind to become so over-ionized that the fast stream 
produced by disk dominated spectra gives way to a 'failed 
wind' structure (Proga 2005). In this case, the material rises 
above the disk duo to UV absorption and starts to accel- 
erate, but as it is exposed to the X-ray source it becomes 
over-ionized, thus preventing continued acceleration, and the 
material eventually falls back down to the disk rather than 
escaping as a wind. This situation can result in quite high 
densities at large scale heights above the disk, complete with 
a complex, circulating velocity field and, crucially, very little 
mass- loss (Proga 2005). Since this occurs directly above the 
inner disk, the velocity dispersion is expected to be large and 
it is possible that magnetic fields may help to stir the mate- 
rial up to even greater scale heights, increasing the covering 
fraction of this material. 



6 CONCLUSION 

GD04 proposed a model for the soft excess, observed in 
the X-ray spectra of many type-1 AGN, based on heavily 
smeared, ionized, absorption. They suggested a tentative 
origin for this absorption in the innermost, accelerating, re- 
gion of an accretion disk wind. We expand on this model 
by including the emission that must be associated with this 
absorption, focusing on the impact that the emission has 
on the total spectrum. Initially, we examine the properties 
of the emission, the absorption and the total spectrum as 
we vary the free parameters in the model, before testing the 
model against data for PG1211+143, a bright nearby quasar 
known to have one of the strongest soft X-ray excesses. 

The main conclusion from this work is that the emission 
associated with heavily smeared, ionized, absorption actu- 



ally helps this model reproduce the soft excess, by fiUing-in 
some of the absorption trough. This removes some of the fine 
tuning inherent in the original pure absorption model and, 
rather than reducing the curvature imposed on the spectrum 
by the model, the emission instead allows the model to re- 
produce the soft excess curvature between 0.6-1.5 keV from 
a wider range of column densities and ionization states, with 
lower velocity smearing. 

We show this specifically by fitting the total smeared ab- 
sorption/emission to the XMM-Newton EPIC PN spectrum 
of PG1211-I-143, as part of a more complex model that ac- 
counts for complex narrow absorption systems and Compton 
reflection from an ionized accretion disk. This model results 
in a good fit to the data (Xre<i~l-08) demonstrating that it 
is capable of reproducing even the strongest soft X-ray ex- 
cesses observed. Furthermore, in comparison with the details 
discussed in GD04, the addition of the emission reduces the 
column density required to fit these data by a factor ^2 and 
reduces the smearing velocity from ~0.28c to ~0.2c. The 
accretion disk reflection component does contribute to the 
spectrum, but it does not produce a strong contribution to 
the soft excess. The reflection is smeared, but not dramati- 
cally so. This is in sharp contrast to the reflection dominated 
models for the soft excesses, which require extreme spin of 
the black hole. 

However, there are significant problems with interpret- 
ing the absorbing material as an accretion disk wind. The 
interpretation as a UV line-driven disc wind, suggested in 
GD04, is especially problematic because such winds tend to 
be produced within ~25° of the equatorial plane (PK04). 
Thus most type-1 AGN are expected to have lines of sight 
which would not intercept much of this material. This di- 
rectly conflicts with the observations that most PG Quasars 
(type-1 AGN) show a soft excess. There is also a more 
generic problem in identifying this material as a wind. The 
required large velocity smearing implies a mass-loss rate 
of Moutfiom ~190 M0 yr^"', much larger than the accre- 
tion rate required to power the observed luminosity of ~1- 
2 M© yr"^ (see also Chevalher et al. 2005). 

We suggest an origin for this material as a 'failed wind', 
i.e., that the central X-ray source is strong enough to over- 
ionize the wind, removing the acceleration through line ab- 
sorption before the material reaches escape velocity. Thus 
the material is initially accelerated, but it cannot escape 
the system and eventually falls back to the disk. The un- 
reasonably large mass-loss rate is then easily circumvented, 
since most of the material simply circulates, falling back and 
rejoining the disk at somewhat larger radii. Modelling such 
material is undoubtedly much more complex than described 
by our very simplistic assumptions on the density and veloc- 
ity structure, but despite this we get a very good fit to the 
shape of the soft excess in PG1211-I-143 from a combination 
of absorption and emission in this material. Better models 
fit to a wide range of type-1 AGN will allow future studies 
to determine if this interpretation is physically plausible. 
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APPENDIX A: DETAILS OF THE EMITTED 
LUMINOSITY 

The equations used for photoionization codes {e.g. 
CLOUDY, XSTAR, TITAN, etc) are complex and are pre- 
sented in detail elsewhere (most notably in the manuals for 
the specific codes). However it is instructive to understand 
some of the dependencies of bound-bound (lines), bound- 
free (RRC) and free-free (ff) interactions in a photoionized 
gas. 

Equations Al, A2 & A4 arc expanded versions of equa- 
tion 2 (Section 2), showing some the complexities of the 
F„-functions, We note that even these equations have been 
considerably simplified, for example we ignore coUisional 
ionization and three body interactions. Each of these equa- 
tions basically contains a sum of the energy released from 
the probable interactions for each ion species, i, in the gas, 
weighted by the fraction of the gas made up by this ion, fi 
(Rybicki & Lightman 1979, Osterbrock 1989). 



Al Radiative Recombination 

The luminosity emitted by radiative recombination transi- 
tions is given by, 



Lrrc = nlV fi I yTtflflc.i.j! / 

i \ X Li./, 



hvi (x) dvi 



(Al) 



Again, is the electron density of the material, fi is the 
fraction of the gas in the specific ion species, i. aRRc.i.x is 
the radiative recombination coefficient for the recombina- 
tion of a free electron to a specific energy level, x, in the 
ion species, i. The integral is the total energy released by 
recombinations into energy level, x, for ion, i, by the total 
distribution of free electrons, where Vmin is the minimum 
frequency of a photon released by the recombination of a 
barely-free electron to energy level, x. 
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A2 Line emission 

Since line emission in a pliotoionized gas relies (to first order) 
on a previous radiative recombination in order for a line 
emitting cascade to occur, the luminosity emitted in line 
transitions is given by, 




rie is the electron density of the material, fi is the fraction 

of the gas in the specific ion species, i, and aRR.c,i is the ra- 
diative recombination coefficient for ion species, i. Pi,t is the 
probability of a particular set of atomic energy level transi- 
tions, t, occurring for the specific ion species, i, and hi>i,t is 
the total energy of the photons released by that particular 
set of atomic level transitions. Here aRRc,i is given by, 

CtRRC,i = <XRRC,i,x (A3) 

X 

We note that the nl dependence of this component is forced 
on it by the requirement that a 2-body (electron-ion) recom- 
bination event occurs prior to the line emission. 

A3 Free-free emission 

The luminosity emitted by free-free integrations is given by, 
Lff = nlvY^fi aff,i (^j hvi dv^ (A4) 

Tie is the electron density of the material, fi is the fraction 
of the gas in the specific ion species, i. ctff^i is the free- 
free coefficient for ion species i and the integral is the total 
energy released by the free-free integrations of ion, i, with 
the total distribution of free electrons. 



material, however this condition then applies throughout the 
shell of material, not just to the surface layers. 



APPENDIX B: PHOTON LOSS THROUGH THE 
SURFACES OF A THICK WIND 

Consider a simple 'toy' model of an accretion disk wind £is 
uniform outflowing gas in a conic section torus with a with 
an opening angle, 0, and a well-defined upper, lower and 
outer boundary. Photons entering the material are repro- 
cessed through absorption, re-emission and scattering and, 
to first order, reprocessed photons axe distributed isotrop- 
ically with each interaction. Photons produced by interac- 
tions near the upper/lower surface of the wind arc capable of 
escaping from the wind through these surfaces, significantly 
changing the amount of energy available for following inter- 
actions in the remaining material (King & Pounds 2003). 

If the wind is radially thin this effect is negligible, 
given the small surface area available for photons to be lost 
through, and so a thin shell (no upper/lower surfaces) is a 
reasonable approximation of a thin wind. This is the case 
presented in the model here. This is not the case for a thick 
wind which, with its considerably larger surface area, will 
not be well approximated by a thick shell of gas. A thick 
wind can be crudely approximated with XSTAR by allow- 
ing the local covering fraction to be, cfracKl. This allows 
photons to escape within each calculational sub-layer in the 
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